Gaussian quantum discord 
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We extend the quantum discord to continuous variable systems and evaluate Gaussian quantum 
discord C'(q) for bipartite Gaussian states. In particular, for squeezed thermal states (STS), we 
explicitly maximize the extractable information over Gaussian measurements: C(q) is minimized 
by a generalized measurement rather than a projective one. Almost all STS have nonzero Gaussian 
discord: they may be either separable or entangled if the discord is below the threshold C(q) = 1, 
whereas they are all entangled above the threshold. We elucidate the general role of state parameters 
in determining the discord and discuss its evolution in noisy channels. 
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Quantum correlations have been the subject of inten- 
sive studies in the last two decades, mainly due to the 
general belief that they are a fundamental resource for 
quantum information processing tasks. The first rigorous 
attempt to address the classification of quantum correla- 
tion from has been put forward by Werner [l[ , who put on 
firm basis the elusive concept of quantum entanglement. 
A state of a bipartite system is called entangled if it can- 
not be written as follows: qab — J^PkQAk ® QBk, where 
QAk and QBk are generic density matrices describing the 
states of the two subsystems. The definition above has an 
immediate operational interpretation: separable states 
can be prepared by local operations and classical com- 
munication between the two parties, whereas entangled 
states cannot. One might have thought that such clas- 
sical information exchange could not bring any quantum 
character to the correlations in the state. In this sense 
separability has often been regarded as a synonymous 
of classicality. However, it has been shown that this is 
not the case @, 0] . A measure of correlations -quantum 
discord- has been defined as the mismatch between two 
quantum analogues of classically equivalent expression of 
the mutual information. For pure entangled states quan- 
tum discord coincides with the entropy of entanglement. 
However, quantum discord can be different from zero also 
for some (mixed) separable state. In other words, clas- 
sical communication can give rise to quantum correla- 
tions. This can be understood by considering that the 
states QAk and QBk above may be physically non distin- 
guishable, i.e. non-orthogonal and thus not all the infor- 
mation about them can then be locally retrieved. This 
phenomenon has no classical counterpart, thus account- 
ing for the quantumness of the correlations in separable 
state with positive discord. Quantum discord has been 
shown to be a property hold by almost all quantum states 
4] and recently attracted considerable attention [jl-fioj]. 
In particular, the vanishing of quantum discord between 
two systems has been shown to be a requirement for the 
complete positivity of the reduced subsystem dynamics 
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While the discord is a fundamental notion allowing for 
the description of the quantumness of the correlations 
present in the state of a quantum system, its evalua- 
tion requires an optimization procedure over the set of 
all measurements on a given subsystem, and thus attack- 
ing the general case is a formidable task. For this rea- 
son, the original definition of the quantum discord [2[ in- 
volved orthogonal measurements, and its evaluation and 
the study of its properties has mainly been restricted to 
final dimensional systems [l3j |. The purpose of this Let- 
ter is to extend the notion of discord to the domain of 
continuous variable systems. In the following, we focus 
our analysis on bipartite systems that are described by 
two-mode Gaussian states and we explore the concept of 
discord within the domain of generalized Gaussian mea- 
surement, i.e. any measurement that may be achieved 
using passive and active linear optics, homodyne detec- 
tion and auxiliary modes prepared in Gaussian states 
14l |15| . We start our discussion by reviewing the main 
ideas at the basis of the definition of the discord. Let 
us consider two classical random variables A and B with 
joint probability pab (fl, b); the total correlations between 
the two variables are measured by the mutual informa- 
tion. The latter may defined by two equivalent expres- 
sions I(A; B) = H(A) + H(B) - H(A, B) and I(A; B) = 
H(A) - H(A\B) = H(B) - H(B\A) where H(X) = 
— J2 x Px( x ) logpx(a^) is the Shannon entropy of the cor- 
responding probability distribution and the conditional 
entropy is defined in terms of the conditional probability 
PA\s{a\b) as H{A\B) = - Y,abPAB(a, b) \ogp A \ B (a\b). 
The idea of quantum discord grows out of the fact that 
the quantum version of the mutual information of a bi- 
partite state qab may be defined in two nonequivalent 
ways. The first is obtained by the straightforward quanti- 
zation of I( A; B), i.e. I (qab) = S(q a ) + S(q b ) ~ S(qab) 
where S(q) = — Tt[q logg] is the Von-Neumann entropy 
of the state q and Qa(b) — r ^B(A)[QAB\i are the partial 
traces over the two subsystems. On the other hand, the 
quantization of the expression based on conditional en- 
tropy, i.e. the extractable information, involves the con- 
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ditional state of a subsystem after a measurement per- 
formed on the other one and this fact has three relevant 
consequences: i) the symmetry between the two subsys- 
tems is broken; ii) this quantity depends on the choice 
of the measurement; iii) the resulting expression is gen- 
erally different from I(qab)- Let us denote by gAk — 
l/pB(k)Ti B [g AB l(g,n k } with p B (k) = Tr AB[QABl®Hk], 
the conditional state of the system A after having ob- 
served the outcome k from a measurement performed on 
the system B. In turn, {Ilfc}, J^k Hfe = I denotes a prob- 
ability operator-valued measure (POVM) describing a 
generalized measurement. The quantum analogue of the 
mutual information defined via the conditional entropy 
is defined as the upper bound J a = supm fe i S(@a) ~ 
^2kPB{k)S(QAk) taken over all the possible measure- 
ment. Finally, the quantum A-discord is defined in terms 
of the mismatch C{qab) — I(qab) ~ Ja(qab)- Analo- 
gously one is led to define the B-discord through the en- 
tropy of conditional states of system B. In the following 
we show that the extractable information J(qab) for two 
modes gaussian states can be maximized over the class of 
Gaussian measurements, and that the mismatch between 
C(qab) is actually minimized by a POVM rather than 
a projective measurement. As we will see it is enough 
to focus on A-discord since the results for the B-discord 
are recovered by a repameterization of the state: from 
now on we refer to A-discord as the discord of the quan- 
tum state g and omit the indication of the subsystem. 
Recently, a different quantity has been introduced [l~2j ]. 
which is essentially a symmetrized version of the discord. 

We start our analysis by proving a general result: 
quantum discord is invariant under local unitary oper- 
ations, i.e. C(Ua ® UbQabU a ® U B ) = C(qab), Vg 
and any choice of the local unitaries. The proof sim- 
ply follows by first noticing that the mutual informa- 
tion I(qab) is written in terms of two- single-system en- 
tropies and thus it is not changed by the action of local 
unitaries. Furthermore, extractable information rewrites 
as J(g) — S(qa) — Y^kP'si^^i^Ak) wnere the primed 
quantities are evaluated using the transformed POVM 
H' k = U B IlkUB- Since this amounts to a repametrization 
of POVMs, which does not change the superior, invari- 
ance is proved. This result is relevant since it allows us 
to focus our analysis on Gaussian states whose covari- 
ance matrix is in a standard form. Indeed, let us now 
consider bipartite Gaussian states i.e., states that can be 

/ A C \ 

characterized by their covariance matrix cr = ( ^ ^ J . 

By means of local unitaries that preserve the Gaussian 
character of the state, i.e. local symplectic operations, 
<j may be brought to the so-called standard form, i.e. 
A = diag(a, a), B = diag(6, b), C = diag(ci, C2). The 
quantities I\ = det A, I 2 = dot B, I 3 — det C, I4 = det a 
are left unchanged by the transformations, and are thus 
referred to as symplectic invariants. The local invariance 
of the discord has therefore two main consequences. On 



the one hand, C(g) may be written in terms of sym- 
plectic invariants only. On the other hand, it allows 
us to restrict to states with <x already in the standard 
form. In particular, while the derivation we give for 
the Gaussian discord is applicable to the general case, 
for the explicit calculations we will focus on the rel- 
evant subclass of states for which cl = — c2, i.e the 
squeezed-thermal states (STS) g — S(r)v\ ® z^S^r), 
where S(r) = e r ( a1b1 - ab ) [ s the two-mode squeezing op- 
erator and Uj = E/c AjX 1 + V fc )-' £ - 1 |fc)(fc|, j = 1,2 are 
chaotic states with Nj average number of thermal pho- 
tons. Using this parametrization we have a — (N r + i) + 
N 1 (l + N r )+N 2 N r , b = (V r + i)+V 2 (l + V r ) + V 1 V r , 
and ci = -c 2 = (1 + Ni + N 2 ) y / N r (l + N r ) where 
N r = sinh 2 r. 

The definition of the Gaussian quantum discord is 
based on the minimization of the mismatch 1(g) — J(g) 
over single mode generalized Gaussian measurements. A 
first class of such POVMs may be written as [3, [HI 
n x = D{X)g M D\X), JdXUx = t, where X is 
a two dimensional real vector and gu a generic zero 
mean Gaussian state whose covariance matrix is ctm = 

( a a Ii with fixed parameters a, 8 £ R + , 7 € R. 
V7 8 J 

If one performs the measurement described by ITx on, 
say, mode B of a bipartite Gaussian state, then the dis- 
tribution of the outcomes p(X) is a bimodal Gaussian 
with covariance matrix (B + ctm), whereas the condi- 
tional state gx of mode A is a Gaussian state of mean 
X T (B + gm) x G t and covariance matrix given by the 
Schur complement a P = A - C(B + cr M ) _1 C T [IEE3]- 
Quantum discord may be written as C(g) = S(qb) — 
5(g)+inf| nx } J r dX p(X) S(gx) and the general form of 
Gaussian quantum discord is: 

C(g) = h{sfh) - h(d-) - h(d+) + inf h[^p~) (1) 

cm 

where, h[x] = (x + |)log(x + |) — (x — |) log(x — |) 
and d± are the symplectic eigenvalues of g, expressed by 
d% = \ [A ± VA 2 - 4/4 J , A = h + I 2 + 2/ 3 - In deriv- 
ing the expression for C(g), we have used two facts: i) 
the entropy of a Gaussian state depends only on the co- 
variance matrix and ii) the covariance matrix up of the 
conditional state does not depend on the outcome of the 
measurement itself. This facts allows for a simplification 
of the minimization required to obtain the final expres- 
sion general expression of the Gaussian discord. Indeed, 
for the relevant case of STS, and for any choice of N\, 
N2 and N r , the minimum of the mismatch 1(g) — J{g) 
is obtained for a — 8 — 1/2, 7 = i.e. when the co- 
variance matrix of the measurement is the identity. This 
corresponds to the coherent state POVM, i.e. to the 
joint measurement of canonical operators, say position 
and momentum, which may realized on the radiation 
field by means of heterodyne detection [l8| . It turns out 
that the same result is obtained even if we generalize the 
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class of Gaussian measurements to include non covariant 
ones Hz = D(X)qm{Y)D^ (X), where now the vector 
Z = (X,Y) includes the no longer fixed parameters of 
the covariance matrix a^j — um(Y)- Indeed, since the 
integrand in inf^n^}. JdZp(Z) S(gz) is always positive 
we have inf { n z } JdZp(Z) S(g z ) > inf {n z } S{g z ) = 
infy h(y/<jp(Y)) and the above results apply for any 
Z. Upon substituting o~m ~* 1/2, we can now explic- 
itly write the Gaussian discord for the generic bipartite 
STS in terms of of symplectic invariants as: 



C(g) = h{y/I 2 )-h(d-)-h(d + )+h( 



l + 2Jh 



)■ 
(2) 

Upon exchanging 1\ f-> I 2 one can pass from the A- 
discord to the £?-discord. 

We are now ready to start our discussion about the prop- 
erties and the operational meaning of Gaussian quantum 
discord. At first we notice that C(g) ^ as far as N r 7^ 0. 
Given that Gaussian states in standard form are separa- 
ble for N r < N1N2/O- + Nx + N2), this confirms that 
for CV Gaussian states there are separable states with 
nonzero discord. Besides, since JV r / <S> c ^ we 
have that bipartite Gaussian states have always nonzero 
Gaussian discord, except when they are product states. 
The same condition characterizes the class of tomograph- 
ically faithful states for reconstruction of quantum oper- 
ations [19j , and this provides an operational meaning for 
the quantum correlations in separable states with posi- 
tive discord. The behavior of C(g) for small and large 

N r is given by C(g) ^ fx(Nx,N 2 )N r and C(g) ^ 
f2{Nx,N 2 ) + f 3 (Nx ^2) log N r respectively, where fx is 
a decreasing function of -/V2 at any fixed value of Ni and 
f 2 , /3 are decreasing functions of both Nx and N 2 . 
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FIG. 1: (color online) Left: Gaussian discord C(g) for STS at 
separability threshold as a function of the ratio Ni /Nt for in- 
creasing values of the total energy Nt of the states (from bot- 
tom to top N T = 1, 5, 10, 50, 10 2 , 10 3 , 10 s ). Right: symplectic 
eigenvalues of the partial transpose d_ versus C(g) for ran- 
domly generated STS. The red region corresponds to symmet- 
ric states. We also report the separability threshold d_ = ~ 
and the corresponding threshold C(g) = 1 (C(g) = 2 log 2— 1 
symmetric STS). 

We now focus our attention on how C(g) relates with 
other meaningful properties of the states. In Fig. [T^ 
we report C(g) at the separability threshold N r — 



Nx N 2 1 ( 1 + Nx + N 2 ) , as a function of the ratio Nx /N T for 
increasing (from bottom to top) values of Nt = a + b—1, 
which is the total energy of the Gaussian state under 
investigation (Nt = Nx + N 2 + 2NxN 2 at separability 
threshold). The plot suggests two important facts. First, 
Gaussian discord is an increasing function of the total en- 
ergy, and is maximized when most of the thermal photons 
are placed on the unmeasured system, thus maximizing 
the purity of the measured one. Second, the Gaussian 
discord for separable states is always smaller than one. 
The existence of a bound has been confirmed numerically 
by the random generation of a large number of bipartite 
Gaussian states in the standard form: in Fig. [TJd, we re- 
port the smaller symplectic eigenvalue d- of the partially 
transposed state, obtained by replacing I 3 — > —I3 in the 
the formula for <i_, as a function of Gaussian discord. 
Since a Gaussian state is entangled iff d- < i we have 
that for < C(p) < 1 we have either separable or en- 
tangled states, whereas all the states with C(g) > 1 are 
entangled. 

The relation between the discord and the entanglement 
can be further clarified by analyzing the case of symmet- 
ric STS, i.e., Nx = N 2 — N s . Here we focus on the 
behavior of C(g) with respect to global purity of the 
state fi = (1 + 2N S )- 2 and <L = e~ 2;s (l + 2N s )/2. A 
first important observation is that for fixed purity C(g) 
turns out to be a growing function of the entanglement, 
whereas at fixed values of d_ the behavior is more in- 
volved. In Fig. [2] we plot C(g)(fi,d-) at fixed val- 
ues of d-. We can distinguish two different cases. For 
non entangled states (<f_ > 1/2), C(g) decreases with 

and it thus is an increasing function of the total en- 
ergy of the state N T = d- - 1 + (4(L/z) _1 . The lim- 
iting value is thus reached at infinite energy and the 
latter is in general given by C(g)(fi — > 0, d_) = (1 + 
2dl) In [(1 + 2dJ)/d-)] - (1 + (1 +2cL) In 2). Therefore, 
for non entangled symmetric states C(g) < 2 log 2 — 1; 
the latter bound is also reported in Fig. [J and defines 
the limit of the red region corresponding to symmetric 
separable states with non zero discord. 
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FIG. 2: (color online) Gaussian discord C(g) for symmet- 
ric STS (jVi = ./V2 = N a ) as a function of the global pu- 
rity fi of the state p and the smaller symplectic eigenvalue 
of its partially transpose d_. Left: separable states; from 
bottom to top d- — 0.5,0.6,0.7,0.8,0.9; C(g) monotonically 
decreases with pi. Right: entangled states ; from bottom to 
top d- = 0.5, 0.4, 0.3, 0.3, 0.1, 0.06284; G{g) monotonically in- 
creases with fi only for d_ > 0.06284. 

As for the entangled states (d- < 1/2), the behavior of 
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C(g) with p is more complex. For states which are highly 
entangled (d_ < 0.06284) the discord decreases (grows) 
monotonically with p (Nt)- Indeed, J(g) ps h(\/Jx) 
when d_ — > 0, i.e, the extractable information is max- 
imized, and C(g) is maximum for pure states. For inter- 
mediate values of the entanglement, d_ € (0.06284,0.5), 
C(g>) has a non monotonic behavior with (Nt). In 
particular, the maximum discord is reached for /i = 1 
(pure states) only for d- > 0.1282, while its minimum 
is reached for intermediate values of p that depend the 
actual value of cL. The overall non monotonic behav- 
ior of C(q) corresponds to a situation in which, at fixed 
value of entanglement, the quantumness of the state as 
measured by the Gaussian discord, varies depending on 
the total correlations present in the state, and conse- 
quently the ordering of the states with respect to their 
quantumness significantly differs by that given by the en- 
tanglement. We also emphasize that by fixing the value 
of d- one also fixes the value of the teleportation fidelity 
F = (1 + 2d_)~ 1 of coherent states [2fj. This means that 
by varying the global purity of the state p shared by Alice 
and Bob, the same fidelity can be achieved with different 
quantum resources as measured by the Gaussian discord. 

We finally address the fundamental issue of the evo- 
lution of quantum discord in noisy channels. Let us 
consider bipartite Gaussian states that evolve accord- 
ing to Lindblad Master equation g—\ J2j ^jMjL[a}g + 
r,-(l + Mj)L[a^]g, which describes the Markovian inter- 
action of the two modes with independent thermal reser- 
voirs, Tj and Mj being the damping factor and the av- 
erage number of thermal photons of the two reservoirs 
respectively. The mapping induced by the ME is Gaus- 
sian and the covariance matrix of the evolved state is 
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2 , M 2 + 5 , M 2 + \) is the covari- 
ance matrix of the reservoir, which also describes the sta- 
tionary state of the system. If <x t= o is in standard form its 
parameters evolve as: a' = a e~ Tlt + (1 — e~ Flt )(Mi + |), 



-r 2 t 
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)(M 2 + i), d = ce -5(ri+r2)t ; 



i.e. 

a' > a, b' > b, and d < c. Since C(g) is a decreasing func- 
tion of a and b and an increasing function of c we have 
that Gaussian discord monotonically decreases in noisy 
channels. On the other hand, it has been shown that the 
decrease should be smooth since an arbitrary Markovian 
evolution can never lead to a sudden disappearance of 
discord Q ■ An open question remains the effect of non- 
Markovian dynamics, which has been proved to produce 
oscillations in the dynamics of Gaussian entanglement 
2lL l22j . We also expect Gaussian discord to increase if 
the two parties interact with a common reservoir [23| . 

In conclusion, in this Letter we have extended the no- 
tion of the discord [2J to continuous variable systems and 
discuss its properties. We have defined the Gaussian dis- 
cord C(g) for two modes Gaussian states and we have 
shown the general analytical procedure to derive it. In 



particular, for the relevant subclass of squeezed thermal 
states (STS), we have shown that the extractable infor- 
mation is maximized by a generalize measurement i.e. 
the coherent state POVM corresponding to heterodyne 
detection. Just as the entanglement, C(g) is invariant un- 
der local unitary operations and it is zero only for (ther- 
mal) product states. For separable states C(g) grows 
with the total energy and it is bounded. Numerical ev- 
idences show that in general C(g) < 1 while analytical 
calculations show that for separable symmetric STS the 
bound reduces to C(g) < 2 In 2 — 1. For symmetric STS 
we have also shown that the behavior of C(g) strongly 
depends on the amount of entanglement present in the 
state: it increases with the total purity only when the en- 
tanglement is large, whereas it shows a richer behavior for 
smaller values of entanglement. Our results pave the way 
for the general discussion about the quantum discord in 
continuous variable systems and for its the experimental 
determination with current technology. 
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